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ORBITAL TRANSFER IN MINIMUM TIME 
W. E. BLEICK,* U. S. Naval Postgraduate School 

1 . Introduction. The problem of orbital transfer discussed 
here is that of scheduling the direction p of constant momentum 
thrust of a rocket, which loses mass at a constant rate, so that 
it transfers to an earth satellite orbit, with known elements of 
time, position and velocity, in a minimum time T after launching 
of the rocket. The launching conditions are assumed to be fixed. 
This situation is illustrated in Figure 1 for the case of a cir- 
cular orbit. The sector angle B at which the rocket enters orbit 
will be called the rendezvous angle. To aid the discussion imagin- 
ary physical rendezvous of the rocket and satellite is assumed to 
occur at this angle. The time of rocket launch to achieve actual 
physical rendezvous can be determined, of course, only after 
both of the unknowns T and B have been found. The problem is set 
up as a calculus of variations problem of the Lagrange type, and 
is solved by an iterative process in which an initial approxima- 
tion to the angle B is estimated. 

A non-rotating Oxy rectangular coordinate system with origin 
at the earth's center is used. The coordinates and velocity com- 
ponents of the rocket and target satellite are denoted by x,y,u,v 
and X,Y,U,V respectively. For simplicity the equations of motion 
of rocket and target will be written in a "non-dimensional" form 
by the use of suitable units. The unit of length is taken as the 

* Supported by the Office of -^aval Research. 
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Fig. 1 . Orbital transfer. 



earth's equatorial radius, R^ = 20,925,000 feet. The unit taken 
for time t is the time required by a hypothetical eart|i satel- 
lite, in equatorial, circular, vacuum, sea level orbit, to 
traverse a sector of one radian. This unit of time is ^R / g = 
13.459 minutes, where g = 32.086 ft. /sec. is the acceleration 
of gravity at the equator. The unit of velocity is then the speed 
of this hypothetical satellite. These units of length and time 
will always be understood, unless other, more conventional units 
are specifically mentioned. 

2. Statement of the problem. The equations of motion of the 
rocket, in terms of the specified units of length and time, are 

(1) ‘ CP 1 = u - g 1 - a cos p = 0 

= v - g^ - a sin p = 0 
cp^ = x — u = 0 
= y - v = 0 

where g 1 = -x/r , g 2 = -y/r , r = x + y , a = cm/gO - mt) , 
where m is the constant fraction of initial gross rocket mass 
lost per unit of time, c is the constant speed of the emitted 
rocket gases, and g is the acceleration of gravity at the equator. 
The fixed initial conditions of the rocket trajectory are taken as 

(2) x(0) = 0, u(0) = V 1 = Vq cos 9 

y(0) = 1 , v( 0) = V 2 = V Q sin 6 . 
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The terminal point of the rocket trajectory is variable with 



(3) x(T) =X(T), y(T) = Y(T) , u(T) = U(T) , v(T) = V(T) . 



It is also assumed that the rocket thrust is turned off abruptly 
at time T. This discontinuity will lead to a trivial steering 
corner in the calculus of variations problem. Note that (1) and 
(3) imply that U(T) = -[X/R 3 ] t = g., (T) and Y(T) = -[y/R 3 ] t = g 2 (T) 
where R = |? X + ^ Y| . 

The problem is to choose the control variable p to effect 



r 

Jo 



orbital transfer with L dt minimized, and to determine -the 
corresponding rocket trajectory. This problem is equivalent to 
the Lagrange calculus of variations problem of requiring the 
integral 



(4) 



f T 

= Jo 



(1 + Xcp-j + ucp 2 + trcp 3 + pcp 4 )dt 



to be stationary. In (4) the equations (l) are regarded as con- 
straints with X(t), li(t) , rr(t) , p(t) introduced as continua of 
Lagrangian multipliers [ 1 ] . If the time coordinate of the varied 
terminal point is taken as T + AT , the vanishing first variation.. 
[2] of I is 

(5) 61 = /o T [X(6u - g lx; 6x - g-| y 6y + a sin p6p) + tt(6x - 6u) 

+ ii(6v - g2 x ^ x - &2y^ ~ a cos + P( 6 Y “ 6v) ]dt 

f T+AT 

+ Jm [l + X6u + n6v - (X cos p + pi sin p)6a]dt = 0 
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where the finite variation 6a, 6u, 6v terms in the integral from 
T to T + AT , created by thrust termination at time T on the un- 
varied trajectory and at T + AT on the varied trajectory, cancel 
On integrating by parts one obtains 

f T . 

(6) 61 = [ \6u + |i6v + tt6x + pdy]^ - J ^ [(X + tt)6u + (|i + p)6v 

+ ( n + g !x X + g 2x^ 6x + (P + Si y X + g 2y^ 6y 
+ a(|-i cos p - X sin p)6p]dt + AT = 0. 

The variations of the dependent coordinates at the variable 

terminal point must be taken [.2] as 

(7) 6u(T) = (U - u) t AT = -(a cos p) T AT , 6x(T) = (X - x) T AT = 0 

6v(T) = (V - v) t AT = -(a sin p) T AT, 6y(T) = (T - y) T AT = 0. 

Substitution of (7) into (6), and application of the fundamental 
lemma of the calculus of variations and the theory of ordinary 
extrema, gives the Euler equations 

(8) X + tt = 0 
£ + P = 0 

" + 8i x X + S 2 x u = 0 
P + g, y * + g 2y U = 0 
tan p = |i/X 

and the transversality condition 

[a(X cos p + n sin p) ] T = 1. 
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The first four homogeneous equations of the Euler equations (8) 
constitute the adjoint system [3] of the system of variation 
equations 

(10) 6cp 1 = 6cp 2 = 6cp 3 = 6cp 4 = 0 

which are the coefficients of X, (i, tt , p in (5) equated to zero. 

The adjoint system has a matrix of coefficients which is the neg- 
ative transpose of that of (10). The last of the Euler equations 
(8) requires that the control variable p be adjusted so that 
a = a(? cos p + "j sin p) , which is proportional to the rocket 
thrust, is continually parallel to the adjoint vector A = i X + j |i. 
The transversality condition (9), which may be written 
(a*A)^ = 1 > 0, requires that a and A have the same sense. Since 
it is only the ratio of pi to X which determines p, it is a trivial 
matter to scale them to satisfy the magnitude requirement of (9). 

A solution of the problem obtained from (1) and (8) guarantees 
a stationary time of transfer. The nature of the problem is such 
that this stationary time is a minimum time. 

3. Numerical solution. There is a constructive aspect of a 
modification of equation (6), first used by Bliss [4] in his work'" 
on differential corrections in ballistics, and applied recently 
by Faulkner [ 5 ] in an iterative fashion in optimum control prob- 
lems. To find the desired modification of (6), assume that a 
solution of the systems (l) and (8) has been obtained, which does 
not necessarily satisfy the terminal conditions (3). Using this 
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solution and holding T constant, consider the variation of the 
vanishing integral 



f T 

(11) Jo ( Xcp 1 + |Jcp 2 + n(:p 3 + pcp 4^ dt = 0 



with the terminal constraints (3) removed, so that the terminal 
variations 6u(T) , 6v(T) , 6x(T) and 6(T) become free. Since X, pi, 
tt, p satisfy the adjoint system, and since there is now no steer- 
ing corner due to thrust termination, one obtains 

T 



(12) [X6u + pSv + tt6x + 



pSy] T = J 0 

-Jo 



.(-X sin p + u cos p)6pdt 



A* ( da/ dp) 6pdt 



where a = a(T cos p + sin p) . Equation (12), which is the 
desired modification of (6), is called the fundamental formula 
by Bliss [ 4 ], but is also known under the generic name of Green's 
formula [3], By the use of (12) it is possible to generate the 
control parameters of a varied trajectory which, hopefully, 
comes closer to satisfying the desired terminal conditions (3). 

To do this, assume that the adjoint system has been solved to 
obtain a fundamental set of four linearly independent solutions 
given by the rows of 

(13) Transpose of B(t) = [^(t) ^(t) n^t) P i ( t ) ] i = 1,2, 3, 4 

where B(0) = I is the identity matrix. The solution A = "i X + j” p. 
of the adjoint system, required to satisfy the last of the Euler 
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equations (8) , is taken as the linear combination 



( 14 ) X = X^ + 1^2 + mX^ 4- nX^ 

U = U 1 + 1U 2 + mU 3 + nU 4 
so that the control angle p is determined by 

(15) tan p = (u 1 + lu 2 + raliq + nu 4 )/(X 1 + 1X 2 + mX^ + nX 4 ) 
and its variation by 

(16) 6p = [ (Xu 2 -tiX 2 )6l + (Xn 3 -|jX 3 )6m + ( Xu 4 -|iX 4 ) 6n]/( X 2 + u 2 ) . 

When (13) and (16) are substituted into (12) there results the 
system of Green's formulae 

(17) [6u 6v 6x 6y] T B(T) = [0 61 6m 6n]A 

where the elements of the matrix A are 

08) a.. = JJ a(Xu. - |iX.) (Xu.. - t iX j )dt/(X 2 + U 2 ) 3/2 . 

The coordinates of the terminal point of the varied trajectory 
at time T + AT may be taken as [u+Au v+Av x+Ax y+Ay]^ where 

(19) [Au Av Ax Ay] T = [6u 6v 6x 6y] T + [u v x y] T AT. 

In an effort to make this new terminal point come closer to satis- 
fying the terminal conditions (3), one may take 

(20) [Au Av Ax Ay] T = [U-u V-v X-x Y-y] T + [U V X Y] t AT. 



- 7 - 



* 



Substitute (19) and (20) into (17) to obtain 

(21) [U-u V-v X-x Y-y] T AT + [0 61 6m 6n]AB -1 (T) = [U-u V-v X-x Y-y] T 

as the system of equations for the determination of AT, 6l, 6m, 6n 
on the varied trajectory. The Faulkner [5] scheme for the numerical 
solution of optimum control problems may nov be stated: Make an 
initial guess for the values of T, 1, m, n; carry out a simultan- 
eous numerical integration of the systems ( 1 ) , (8) and (18) using 

the control variable of (15); solve the system (21) for AT and 

% 

the changes in the control parameters; iterate until convergence 
is obtained. This program may be carried out in a matter of seconds 
on modern digital computers. 

To give an example of this control optimization in the present 
problem a circular satellite orbit was assumed of radius R = 

| i X + Y | = 1 .075699, corresponding to an altitude of 300 statute 
miles above sea level. For this orbit |T U + j* V| = 1 /YR and 
| i U + j V | = 1/R . The assumed rocket launching velocity was 
Vq = |T V 1 + j* | = 0.585402, launch angle 0 = 0.928084 and ren- 
dezvous sector B = 0.153840 as in Figure 1. Also assumed were 
c = 10000.9 ft. /sec. and m = 0.00360583 sec. ^ The odd appearance 
of these figures is related to the difficulty, explained in the 
next section, of obtaining the initial "guesses" T = 0.289725, 

1 = -0.223125, m = -29.9875 and n = 19.0847 . With this input the 



process converged in four iterations to the seven significant 
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Pig. 2. Trajectory and thrust directions. 



figure results T = 0.2894592, 1 = 0.1840054, m = -108.94383, 
n = 67.95886 and B = 0.1536015 . Figure 2 shows the resulting 
trajectory and rocket thrust directions for equal time intervals, 
excepting the interval terminating in transfer, which is reduced 
by one half. 

4. The initial guesses. The domain of convergence of the 

iteration scheme of the last section appears to be rather limited 

in the present control problem, requiring initial guesses for 

T, 1, m, n which make [U-u V-v X-x Y-y]^, small in (21). In 

* 

the numerical example given here, where the angle B is small and 
radius R nearly unity, the equations (l) can be linearized and an 

exact solution of the linearized transfer problem used to supply 

/ 

the input guesses required to solve the non-linear problem. The 
gravitational terras g^ and g 2 of (l) were replaced by the linear 
terms in their Tfaylor expansions at x = 0, y = 1. The systems (1) 
and (8) then become uncoupled. Their solutions, using (2), are 



( 22 ) 





cos p sin (t-w)dw + sin t 

cos p cos (t-w)dw + cos t 

sin p sinh "f2(t-w)dw + sinh f2t - (l/Y^)cosh f2t + 3/V2 
sin p cosh Y^(t-w)dw + V 2 cosh ^2 1 - (l/"f2)sinh T2t 



tan p = [l cosh Y^t - (n/V^2)sinh V2t]/(cos t - m sin t) . 
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It will be of no avail to attempt to solve (22) and (3) for 
T, 1, m, n by the Newton-Raphson method, since the Newton-Raphson 
equations are in fact the system (21) of the poorly convergent 
iterative routine of the last section. A substitute procedure 
of solving (22) and (3) for , V 2 > m and B was used. When the 
transfer problem has been solved for a sufficient number of such 
sets (V.j , , m, B) , a basis will be at hand for obtaining desired 

rocket launching conditions by interpolation or extrapolation.. 

It can now be revealed that the numerical example of the last 

section really had its genesis in the assumptions B = tt/16, 0 = tt/4, 

♦ -1 

Vq = 0.65, m = 0.0025 sec. and c = 10000 ft. /sec. It was 

estimated that T = 0.3, 1 = 1.25, m = -1.0 and n = 9.5 would 

satisfy (22) and (3). Using a cluster of five closely spaced 

points around (T,l,m,n), a single application of regula falsi [6] 

was made in the hope of improving these values of T, 1, m, n. The 

particular cluster chosen gave the output values T = 0.289725, 

1 = -0.223125, m = -29.9875, n = 19.0847, which will be recognized 

as the initial guesses of the last section. The output^ residuals were 

[X-x, Y-y, U-u, V-v] T = [0.041389, -0.020160, 0.212093, -0.221912]. 

These residuals were then processed by the linear system 

f T 

(23) (X-x)^, = AY^ sin T - Y(T)dB + Ac Jq a cos p sin (T-w)dw/c 

(U-u)m = AV. cos T - V(T)dB + Ac f n a cos p cos (T-w)dw/c 

_ J fT 

1/2 (Y-y) ^ = AV 2 sinh V2T + V2X(T)dB + Ac Jq a sin p sinh Y*2(T-w)dw/c 

(V-v) T = AV 2 cosh Y2T + U(T)dB + Ac Jq a sin p cosh ^2(T-w)dw/c, 
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derived from (22) and (3) for the case of a circular orbit, to 
obtain c + Ac, launching conditions + AV^ and "S ' 2 + , and ren 

dezvous sector B + dB which will give small residuals. The value 
of c + Ac turned out to be far from the desired 10000 ft. /sec. 
Instead of changing the value of c, the transformation 

(24) 10000 In [l - (m + Am)T] = (c + Ac) In (1 - mT) , 



(22) and (23). Seven iterations produced the input values of , 



5. The problem of two fixed end points. The problem in which 
the conditions at both ends of the rocket trajectory are fixed, 
or the problem in which rocket and target satellite achieve 
actual physical rendezvous at a specified angle B, can also be 
solved by the methods presented here, but with greater converg- 
ence difficulties. 

Acknowledgement . I am indebted to my colleague, Prof. Frank 
D. Faulkner, for suggesting the problem of this paper and for 
discussions of difficulties encountered. 
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CDC-1604 ASSEMBLY LANGUAGE COMPUTER PROGRAM 
Set content of IC1 00+1 « SLJ 0 L3440 ZRO 0 00000 



KO 



K 1 0 



K20 



K30 
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REM 

REM 

REM 

REM 

SIL 
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1 


NT 


COMPUTE LINEARIZED TRAJECTORY 
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(Q) = DELTA DEPARTURE ANGLE 
STORE NO. OF TRAJS. 


STA 


0 


THETATEM 


STORE 1ST DEPART ANGLE 


STQ 


0 


DELTHETA 


STORE DELTA DEPART ANGLE 


LDA 


0 


UNITY 




FDV 


0 


R 


RECIP. ORBIT RADIUS 


SLJ 


4 


SQROOT 


RET. JUMP TO SQUARE ROOT 


STA 


0 


V 


STORE TARGET SPEED 


FMU 


0 


V 




FOV 


0 


R 




STA 


0 


VSQDR 


STORE TARGET ACCEL. 


LDA 


0 


B 


TARGET SECTOR AT RENDEVOUS 


SLJ 


4 


TRIG+70 


RET. JUMP TO TRIGONOMETRIC SR. 


STA 


0 


COSB 




FMU 


0 


R 




STA 


0 


CAPYFIN 


TARGET COOR. AT RENDEVOUS 


LDA 


0 
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SIU 


4 


K20 + 6 
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SLJ 


4 


TRIG 
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0 
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0 
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STA 


0 


CAPXFIN 
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LDA 


0 
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0 
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STA 


0 
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0 


V 
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0 
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0 
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0 
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0 
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STA 


0 


CAPUDFN 
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0 
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FDV 


0 
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4 


SQROOT 
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FMU 


0 
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STA 


0 
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0 
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0 
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STA 


0 


A 
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STA 


0 


THETA 
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4 


0 
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4 
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0 
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STA 


0 
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4 
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0 
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0 
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SLJ 


0 


0 




LDA 


3 


T 1 




STA 


0 


TAUFIN 




LAC 


0 


OMEGA 




FMU 


0 


TAUFIN 




FAD 


0 


UNITY 




STA 


0 


STOR1 




FDV 


0 


OMEGA 




STA 


0 


STOR2 




FMU 


0 


ST0R2 




FMU 


0 


STOR2 




FMU 


0 


STOR2 




FMU 


0 


ST0R2 




FMU 


0 


ERRCOF 


K 1 30 


FDV 


0 


TAUFIN 




STA 


0 


STOR3 




LDA 


0 


STOR1 




SLJ 


4 


LOG+66 




FMU 


0 


ST0R3 




SLJ 


4 


SQROOT 




SLJ 


4 


SQROOT 




ZRO 


0 


0 




STA 


0 


DELTAU 




SLJ 


4 


EQUAL 




STA 


0 


DEL T 1 




ENA 


0 


4 




STA 


0 


RUNKT 1 




ENI 


1 


14 




ENA 


0 


0 




STA 


1 


TAU 1 


K 1 40 


I JP 


1 


/-I 




ENI 


0 


0 




LDA 


0 


TAUFIN 




SLJ 


4 


TRIG 




S * A 


0 


SINT 




FMU 


0 


VI 




FSB 


0 


CAPXFIN 



END INNER LOOP ( OR SLJ 0 



DEPARTURE THRUST ANGLE 
MOVE EL5, EM5 AND EN5 



(-VT) 

(-VT/R ) 

B-( VT/R) 

S = INITIAL TARGET SECTOR 
ALARM REENTRY FROM LOG S.R. 
DECIMAL DUMP WITH PANEL 
PARAMETER WORD 

DECIMAL DUMP 

PARAMETER WORD 

OCTAL DUMP 

PARAMETER WORD 



SET NEW THETA 

END OUTER LOOP. (SEE K440) 
HALT. START RENDEVOUS PROG. 
ENTER DELXYUV SUBROUTINE 



( STOR1 ) = 1 - WT 

(ST0R2) = 11-WT)/W 

SQUARE 

CUBE 

FOURTH 

FIFTH 

MULT. BY -15.B-12 



1 — WT 

LOG(l-WT) TO BASE E 



COMPUTED TAU INCREMENT 
EQUALIZE DELTAU AND SET B2 
SET EQUALIZED DELT1 

SET GlLL FOR 4- INTEGRALS 



CLEAR RUNKT1 AREA 



SIN T 

VI SINT 
VI SINT - X 



- 15 - 



L3440 ) 



K 1 50 



K 1 60 



K 1 70 



K200 



STA 


0 


DELXF IN 




L DA 


0 


TAUFIN 




SLJ 


4 


TRI G+70 


COS T 


STA 


0 


COST 




FMU 


0 


V 1 


VI COST 


FSB 


0 


CAPUFIN 


VI COST - U 


STA 


0 


DELUFIN 




LOA 


0 


TAUFIN 




FMU 


0 


R00T2 


T TIMES SQ. ROOT 2 


SLJ 


4' 


EXP+107 




ZRO 


O' 


0 




STA 


0 


C0SHTR2 


STORE HYP. COSINE 


FMU 


0 


ROOTHALF 


C0SHTR2/SQ.R00T 2 


STA 


0 


ST0R1 




LDA 


0 


TAUFIN 




FMU 


0 


R00T2 




SLJ 


4 


EXP+71 




STA 


0 


SINHTR2 


STORE HYP. SINE 


FMU 


0 


V 2 


V 2 SINHTR2 


FSB 


0 


ST0R1 




FAD 


0 


THREDR2 


ADD { 3/SQ.ROOT 2) 


STA 


0 


STOR 1 




LAC 


0 


CAPYFIN 


( -Y ) 


FMU 


0 


R00T2 




FAD 


0 


STOR 1 




STA 


0 


DELYFIN 




LDA 


0 


C0SHTR2 




FMU 


0 


V 2 




STA 


0 


STOR 1 


( STOR 1 ) = V2 C0SHTR2 


LAC 


0 


SINHTR2 


( -SINHTR2 ) 


FMU 


0 


ROOTHALF 




FAD 


0 


STOR 1 




FSB 


0 


CAPVFIN 




STA 


0 


DEL VF I N 




LDA 


3 


EL 1 




FMU 


3 


ELI 


ELI SQUARED 


FAD 


0 


UNITY 


1 + ELI SQ. 


SLJ 


4 


SQRGOT 




ZRO 


0 


0 




STA 


0 


HYP 




LDA 


0 


A 




FDV 


0 


HYP 




STA 


0 


STOR 1 


( STOR I ) = A/HYP = A COS P 


FMU 


0 


SINT 




STA 


0 


SINDOT 


A COSP SINT 


LDA 


0 


STOR 1 


A COSP 


FMU 


0 


COST 




STA 


0 


COSDOT 


A COSP COST 


LDA 


0 


STOR 1 


A/HYP 


FMU 


3 


ELI 


A SINP 


FMU 


0 


SINHTR2 




STA 


0 


SINHDOT 


A SINP SINHTR2 


LDA 


0 


STOR 1 


A/HYP 


FMU 


3 


ELI 


A SINP 


FMU 


0 


C0SHTR2 




STA 


0 


COSHDOT 


A SINP C0SHTR2 


SLJ 


4 


ALPHA 


SET UP GILL ROUTINE 


ZRO 


0 


RUNKT 1 


PARAMETER WORD 


ZRO 


0 


DERI V 1 




IJP 


2 


/+1 


JUMP IF B2 NOT = ZERO 


SLJ 


0 


/+3 


B2=0. END OF INTEGRATION 


ENI 


0 


0 


B2 NOT = ZERO 


SLJ 


4 


ALPHA+ l 


INTEGRATE AGAIN 


ENI 


0 


0 




SLJ 


0 


/-2 


JUMP TO TEST B2 


LDA 


0 


DELXFIN 




FAD 


0 


SIN 




STA 


0 


DELXFIN 




LDA 


0 


DELUFIN 




FAD 


0 


COS 





- 16 - 






K210 



DERIV1 



K220 



K230 



K240 



K25 0 



STA 
L DA 
FAD 
FDV 
STA 
LDA 
FAD 
STA 
SLJ 
LAC 
FMU 
FAD 
STA 
LDA 
FDV 
STA 
LDA 
FSB 
STA 
ENI 
SLJ 
STA 
LDA 
SLJ 
ZRO 
STA 
LDA 
FMU 
SLJ 
STA 
LDA 
FMU 
SLJ 
STA 
LDA 
SLJ 
ZRO 
FMU 
STA 
LDA 
SLJ 
FSB 
STA 
FMU 
STA 
LDA 
FMU 
SLJ 
ZRO 
FMU 
FDV 
STA 
LDA 
FMU 
SLJ 
FMU 
FS8 
STA 
FMU 
FAD 
SLJ 
STA 
LDA 
FDV 
STA 
LDA 
FDV 
STA 
FMU 
FMU 



0 


DELUFIN 




0 


DELYFIN 




0 


S I NH 




0 


ROOT2 




0 


DELYFIN 




0 


DEL VF IN 




0 


COSH 




0 


DEL VF IN 




0 


X YUV 


TO EXIT OF SUBROUTINE 


0 


TAU1 




0 


OMEGA 


(-WT) 


0 


UNITY 




0 


UNMINWT 


1-WT 


0 


A 




0 


UNMINWT 




0 


ADOMINWT 


A/( 1-WT) 


0 


TAUFIN 




0 


TAUl 


(STOR 1 ) * TAUFIN - TAU 


0 


STOR1 


0 


0 




4 


TRIG 


(SINT) * SINE OF (STOR 


0 


SINT 


0 


STOR1 




4 


TRIG+70 




0 


0 




0 


COST 




0 


STOR1 




0 


ROOT2 




4 


EXP+71 




0 


SINHTR2 




0 


STOR1 




0 


ROOT 2 




4 


EXP + 1 07 




0 


C0SHTR2 




0 


TAU 1 




4 


TRIG 


SIN TAUl = -LAM3 


0 


0 




3 


EMI 




0 


ST0R2 


( S TOR 2 ) = - EMI LAM3 


0 


TAU 1 




4 


TRIG+70 


COS TAUl = LAM1 


0 


ST0R2 


LAM 1 + EMI LAM3 


0 


LAM 


LAM = LAM 1 + EMI LAM3 


0 


LAM 




0 


ST0R3 


{ ST0R3 ) = LAM SQUARED 


0 


TAUl 


0 


ROOT2 




4 


EXP+7 1 


SINH TAU1R2 


0 


0 




3 


ENI 




0 


ROOT2 


( - ENI MU4 ) 


0 


STOR2 


. (STOR 2) = - ENI MU4 


0 


TAUl 




0 


ROOT2 




4 


EXP+107 


COSH TAU1R2 = MU2 


3 


EL 1 


ELI MU2 


0 


STOR2 


ELI MU2 + ENI MU4 = MU 


0 


MU 




0 


MU 




0 


STOR3 


LAM SQ. + MU SQ. 


4 


SQROOT 


0 


HYP 




0 


MU 




0 


HYP 




0 


ST0R3 


( STOR3) = SINP 


0 

0 


LAM 

HYP 




0 


ST0R4 


( STOR 4 ) « COSP 


0 


SINT 


COSP SINT 


0 


ADOMINWT 



- 17 - 





STA 


0 


SINDOT 






LDA 


0 


ST0R4 


COSP 




FMU 


0 


COST 


COSP COST 




FMU 


Q 


ADQMINWT 






STA 


0 


COSDOT 






LDA 


0 


ST0R3 


SINP 




FMU 


0 


SINHTR2 


SINP SINHTR2 




FMU 


0 


ADO M I N WT 






STA 


0 


SINHDOT 






LDA 


0 


STOR3 


SINP 




FMU 


0 


C0SHTR2 


SINP COSHTR2 




FMU 


0 


ADOMINWT 






STA 


0 


COSHOOT 


EXIT FROM DERIV1 


K260 


SLJ 


0 


ALPHA+2 




ZRO 


0 


0 




EQUAL 


SLJ 


0 


0 


S.R. TO EQUALIZE DELTAU 


LDA 


0 


TAUFIN 






FDV 


0 


DELTAU 


CONVERT DELTAU TO FIXED PT. 




SLJ 


4 


FIXIT 




STA 


0 


INTSIGN 


STORE INTEGRAL PART 




STQ 


0 


FRACSIGN 


STORE FRACTIONAL PART 




AJP 


1 


/ + 5 


JUMP IF A NOT = ZERO 




QJP 


0 


/+3 


A=0. JUMP IF Q = 0 




LDA 


0 


TAUFIN 


A=0 . Q NOT = 0 




STA 


0 


DELTFLSG 






ENI 


2 


1 






SLJ 


0 


/+ 1 2 






ENA 


0 


0 


A = Q = 0 




STA 


0 


DELTFLSG 




K27 0 


ENI 


2 


0 






SLJ 


0 


/+10 






AJP 


2 


/+1 


A NOT 0. JUMP IF A GREATER THAN 0 




SCM 


0 


MASK 


A LESS THAN 0. ABSO INTEGER TO A 




INA 


0 


1 


ABSO. INTEGER + 1 TO A 




SAU 


0 


/+] 






ENI 


2 


0 


ABSO. INTEGER + 1 TO 82 




SCA 


1 


2057 






LRS 


0 


57 






ENA 


1 


0 






LLS 


0 


44 






STA 


0 


FLOIPLON 






LDA 


0 


TAUFIN 






FDV 


0 


FLOIPLON 






STA 


0 


DELTFLSG 






ENI 


0 


0 




K300 


SIL 


2 


FIXIPLON 


USED IN FINAL PROG. 




LDA 


0 


DELTFLSG 






SLJ 


0 


EQUAL 


TO EXIT OF S.R. 




ZRO 


0 


0 




FIXIT 


BSS 




44 


B 1 = 5 B6 = 65302 P=460 


R 


DEC 




1.075698925 NON-DIMEN. ORBIT RADIUS 


V 


BSS 




1 


NON-DIMEN. TARGET SPEED 


VSQDR 


BSS 




1 


NON-DIMEN. TARGET ACCEL- 


B 


OCT 




1775622077325042 TARGET SECTOR AT RENDEVOUS = 


COSB 


BSS 




1 




SINB 


BSS 




1 




REARTH 


DEC 




20925000. 


EARTH EQUATOR RAD. IN FT. 


GACCEL 


DEC 




32.086 


EQUATOR GRAVITY IN FT/SEC/SEC 


RMS 


BSS 




1 




OMEGA 


BSS 




1 


NON-D I'ME N. MASS LOSS RATE 


THETATEM 


BSS 




1 




VI 


BSS 




1 


NON-DIMEN. INIT. HORIZONTAL VEL. 


V2 


BSS 




1 


NON-DIMEN. INIT. VERTICAL VEL. 


ERRCOF 


DEC 




- 15. B- 12 


(-15/4096 DEC.) USEO TO FIND DELTAU 


ROOTHALF 


OCT 




2000552023631500 


SINT 


BSS 




1 




COST 


BSS 




1 




SINHTR2 


BSS 




1 




COSHTR2 


BSS 




1 




UNMINWT 


BSS 




1 


STORE FOR ( 1-WT ) 



18 



EL 1 INIT 

EM UNIT 

6N1 INIT 

T 1 INIT 

EL2INIT 

EM2INI T 

EN2INIT 

T2INIT 

EL3 INIT 

EM3INIT 

EN3INIT 

T3INIT 

EL4 INIT 

EM4INIT 

EN4INIT 

T4INIT 

ELS INIT 

EM5INIT 

EN5INIT 

T5INIT 

RUNKT 1 

DELT 1 

TAU 1 

SINDOT 

SIN 

COS DOT 
COS 

SINHOOT 

SINH 

COSHDOT 

COSH 

K440 



EXP 

SQROOT 

TRIG 

ADCMINWT 

R0GT2 

TWO 

THREE 

MASK 

TENGRAND 

COMPRMS 

8EGIN 



LIO 



BSS 




2 








OCT 




200 1500400000000 


DEC. VALUE = 


+ 1 .25 1953125 


OCT 




5776377777777777 


DEC. VALUE = 


- 1.000000000 


OCT 




2004460000000000 


DEC. VALUE = 


+9.500000000 


OCT 




17 76463 146314631 


OEC. VALUE = 


+0.3000000000 


OCT 




200 150030000001 0 


DEC. VALUE = 


+1.251464844 


OCT 




5776377777777777 


DEC. VALUE = 


-1.000000000 


OCT 




2004460100000000 


DEC. VALUE = 


+9.503906250 


OCT 




1776461700000000 


DEC. VALUE = 


+0.2987060547 


OCT 




200 1500200000000 


DEC. VALUE = 


+1.250976562 


OCT 




5776377777777777 


OEC. VALUE = 


-1.000000000 


OCT 




2004460200000000 


OEC. VALUE = 


+9.507812500 


OCT 




1776462000000000 


DEC. VALUE = 


+0.2988281250 


OCT 




2001500100000000 


DEC. VALUE = 


+1.250488281 


OCT 




5776377707007777 


DEC. VALUE = 


-1.000434756 


OCT 




2004460300000000 


DEC. VALUE = 


+9.51 1718750 


OCT 




1776464300000000 


DEC. VALUE = 


+0.3011474609 


OCT 




200 1500000000000 


DEC. VALUE = 


+1.250000000 


OCT 




5776377777077077 


OEC. VALUE = 


-1.000006689 


OCT 




2004460400000000 


DEC. VALUE = 


+9.515625000 


OCT 




1776465000000000 


DEC. VALUE = 


+0.3017578125 


BSS 




1 








BSS 




1 








BSS 




1 








BSS 




1 








BSS 




2 








BSS 




1 








BSS 




2 








BSS 




1 








BSS 




2 








BSS 




1 








BSS 




3 








ENI 


3 


23 


RESET 


LMNT INIT. (SEE K120) 


ENI 


0 


0 








LDA 


3 


EL 1 








STA 


3 


ELI INIT 








IJP 


3 


/-I 








LDA 


0 


NT 








A JP 


1 


K20+5 








SLS 


0 


BEGIN 


HALT. 


START RENDEVOUS PROG. 


BSS 




3 








BSS 




1 17 


B 1 = 3 


B6=65447 P=460 


BSS 




65 


B 1=6 


86=65566 P=460 


BSS 




1 16 


B 1 =7 


86=65653 P=460 


BSS 




1 


STORE 


FOR A/ ( 1-WT) 




OCT 




2001552023631500 






OCT 




2002400000000000 






OCT 




2002600000000000 






OCT 




7777777777777777 






DEC 




lOOCO. 








OCT 




1754400000000000 






REM 




12 


RENDEVOUS PROG. STARTS HERE 


ENI 


6 


SET NUMBER OF ITERATIONS 


ENI 


0 


0 








ENI 


0 


0 








ENI 


0 


0 








LDA 


0 


V 








FDV 


0 


R 


V/R = 


W 




STA 


0 


W 


STORE 


TARGET ANGULAR VEL. 


ENI 


0 


0 








ENI 


0 


0 








ENI 


0 


0 








ENI 


0 


0 








ENA 


0 


35 




« 




STA 

ENI 

ENI 


0 

0 

0 


RUNKT 

0 

0 


SET GILL S.R. FOR 29 DIFF.EQUS. 








ENI 


1 


0 








ENA 


0 


0 


CLEAR 


RUNKT+17 AREA 




STA 


1 


RUNKT +17 







- 19 - 



L20 



L30 

RECUR 



LAMMU 



L50 



STA 


1 


RUNKT+20 






INI 


1 


2 






ISK 


1 


1 12 






SLJ 


0 


/-?. 






LDA 


0 


UNITY 




L AM 1 { 0 ) = 1 


STA 


0 


RUNKT+20 


SET 


STA 


0 


RUNKT+37 


SET 


MU2 (0>=1 


STA 


0 


RUNKT+56 


SET 


P I 3 ( 0 ) = 1 


STA 


0 


RUNKT+75 


SET 


RH04 ( 0 )= 1 


STA 


0 


RUNKT+25 


SET 


P 1 1 DOT ( 0 )= 1 


LAC 


0 


UNITY 






STA 


0 


RUNKT+47 


SET 


L AM3D0T ( 0 )= 


STA 


0 


RUNKT+66 


SET 


MU4DOT(0)=- 


LAC 


0 


TWO 






STA 


0 


RUNKT+44 


SET 


RHO2D0T ( 0 ) = 


ENI 


0 


0 






ENI 


0 


0 






ENI 


0 


0 






ENI 


0 


0 






ENA 


0 


0 






STA 


0 


RUNKT+2 


SET 


TAU = 0 


STA 


0 


RUNKT+4 


SET 


X ( 0 ) =0 


LDA 


0 


VI 






STA 


0 


RUNKT+3 


SET 


XDOT ( 0) = V 1 


STA 


0 


RUNKT+12 


SET 


U( 0)=V1 


LDA 


0 


V 2 






STA 


0 


RUNKT+6 


SET 


YDOT { 0 ) =V2 


STA 


0 


RUNKT+ 15 


SET 


V( 0)=V2 


LDA 


0 


UNITY 






STA 


0 


RUNKT+7 


SET 


-< 

o 

It 


ENI 


0 


0 






SLJ 


0 


/+5 


BYPASS 


SIU 


6 


/+1 






SLS 


1 


/ + 1 


STOP SWITCH 


ENI 


6 


0 






ENI 


1 


147 






LDA 


1 


SAVE 






STA 


1 


RUNKT 






IJP 


1 


/-I 






ENI 


0 


0 






SLJ 


4 


EQUAL 






ZRO 


0 


0 






STA 


0 


RUNKT+1 


STORE EQUALIZED 


SLJ 


4 


LAMMU 


RET. 


■ JUMP TO OBT 


SLJ 


0 


L 70 


BYPASS LAMMU 


ZRO 


0 


0 






SLJ 


0 


0 


ENTRANCE TO S.R 


SIL 


1 


L 60+3 


SAVE 


: B 1 


SIU 


3 


L 60+4 


SAVE 83 


LDA 


0 


RUNKT+20 


LOAD LAM 1 


STA 


0 


ST0R3 


L AMI 


TO ST0R3 


ENI 


0 


0 






ENI 


1 


0 


SET 


B 1=0 


ENI 


3 


0 


SET 


B 3=0 


LDA 


1 


EL 






FMU 


3 


RUNKT+34 


ELILAM2) 


FAD 


0 


ST0R3 






STA 


0 


ST0R3 






INI 


3 


14 






ENI 


0 


0 




' 


ISK 


1 


2 






SLJ 


0 


/-i 






STA 


0 


LAM 


STORE LAMBDA 


FMU 


0 


LAM 


LAM 


SQUARED 


STA 


0 


ST0R1 


LAM 


SQUARE TO S 


LDA 


0 


RUNKT+23 


LOAD MU1 


STA 


0 


ST0R3 






ENI 


3 


0 






LDA 


1 


EL 






FMU 


3 


RUNKT+37 


EUMU2) 



- 20 - 



L60 



L70 



DERIV 



L 1 1 0 



L 120 



FAD 


0 


ST0R3 




STA 


0 


ST0R3 




INI 


3 


14 




ENI 


0 


0 




ISK 


1 


2 




SLJ 


0 


/- 3 




STA 


0 


MU 


STORE MU 


FMU 


0 


MU 


MU SQUARED 


FAD 


0 


ST0R1 


LAM SQ + MU SQ 


STA 


0 


HYP 2 


HYP2 = LAM SQ + MU SQ 


SLJ 


4 


SQROOT 




ZRO 


0 


0 




STA 


0 


HYP 


HYP = SQ ROOT OF HYP2 


FMU 


0 


HYP2 




STA 


0 


HYP 3 


HYP3 = HYP2 TIMES HYP 


ENI 


1 


0 


RESTORE B 1 


ENI 


3 


0 


RESTORE B3 


SLJ 

BSS 


0 


LAMMU 

3 


TO EXIT OF S.R. 


L DA 


0 


A 


ENTER FROM L30+7 


FMU 


0 


LAM 




FDV 


0 


HYP 




STA 


0 


RUNKT+11 


SET UDOT(O) = A( LAM) /HYP 


L DA 


0 


A 




FMU 


0 


MU 




FDV 


0 


HYP 




FSB 


0 


UNITY 




STA 


0 


RUNKT+14 


SET VDOT(O) 


ENA 


0 


/ +2 


CURRENT ADORESS+2 TO A 


SAL 


0 


L 220-4 


SET EXIT FROM DERIV. PROG. 


SLJ 


0 


L 150-3 


JUMP TO SET AIJDOT(O) 


ENA 


0 


ALPHA+2 




SAL 


0 


L22C-4 


REPAIR EXIT FROM DERIV. PROG. 


ENI 


0 


0 




SLJ 


0 


L 220-3 


BYPASS DERIV. PROG. 


LAC 


0 


OMEGA 




FMU 


0 


RUNKT+2 


(-WT) 


FAD 


0 


UNITY 


1-WT 


STA 


0 


UNMINWT 


STORE (1-WT) 


LDA 


0 


A 




FDV 


0 


UNMINWT 




STA 


0 


ADOMINWT 


STORE A/M-WT) 


ENI 


0 


0 




LDA 


0 


RUNKT+12 




STA 


0 


RUNKT+3 


SET XDOT = U 


LDA 


0 


RUNKT+15 




STA 


0 


RUNKT+6 


SET YDOT = V 


LDA 


0 


RUNKT+4 


LOAD X 


FMU 


0 


RUNKT+4 


X SQUARED 


STA 


0 


STOR2 


X SQ TO ST0R2 


LDA 


0 


RUNKT+7 


LOAD Y 


FMU 


0 


RUNKT+7 


Y SQUARED 


FAD 


0 


STOR2 


X SO + Y SQ 


STA 


0 


STOR1 


(X SQ + Y SQ) TO ST0R1 


SLJ 


4 


SQROOT 


FMU 


0 


STOR1 




STA 


0 


R32 


STORE RAOIUS CUBED 


FMU 


0 


STOR1 




STA 


0 


R52 


STORE RADIUS FIFTH 


ENI 


1 


0 




SLJ 


4 


LAMMU 


NOTE B 1 SAVED BY LAMMU 


LAC 


I 


RUNKT+4 


(-X) TO A 


FDV 


0 


R32 


(-X/R32) TO ST0R1 


STA 


0 


STOR1 


LDA 


0 


ADOMINWT 


FMU 


1 


LAM 




FDV 


0 


HYP 




FAD 


0 


STOR1 




STA 


1 


RUNKT+ 1 1 


SET UOOT(T) 


INI 


1 


2 



21 



L130 



L 1 40 



L 1 50 



L 1 60 



ENI 


0 


0 


I SK 


1 


5 


SLJ 


0 


L 12C-3 


ENI 


0 


0 


ENI 


3 


0 


LAC 


3 


RUNKT+26 


STA 


3 


RUNKT+ 17 


LAC 


3 


RUNKT+3 1 


STA 


3 


RUNKT+22 


LDA 


0 


RUNKT+7 


FMU 


0 


RUNKT+7 


FSB 


0 


ST0R2 


FSB 


0 


ST0R2 


FMU 


3 


RUNKT+20 


STA 


0 


ST0R1 


LAC 


0 


THREE 


FMU 


0 


RUNKT+4 


FMU 


0 


RUNKT+7 


FMU 


3 


RUNKT+23 


FAD 


0 


ST0R1 


FDV 


0 


R52 


STA 


3 


RUNKT+25 


LAC 


0 


RUNKT+7 


FMU 


0 


RUNKT+7 


FMU 


0 


TWO 


FAD 


0 


ST0R2 


FMU 


3 


RUNKT+23 


STA 


0 


ST0R1 


LAC 


0 


THREE 


FMU 


0 


RUNKT+4 


FMU 


0 


RUNKT+7 


FMU 


3 


RUNKT+20 


FAD 


0 


ST0R1 


FDV 


0 


R52 


STA 


3 


RUNKT+30 


INI 


3 


13 


ENI 


0 


0 


ISK 


3 


57 


SLJ 


0 


L 13C-4 


LAC 


0 


OMEGA 


FMU 


0 


RUNKT+2 


FAD 


0 


UNITY 


STA 


0 


UNMINWT 


LDA 


0 


A 


FDV 


0 


UNM INWT 


STA 


0 


ADOMINWT 


ENI 


0 


0 


LAC 


0 


MU 


FMU 


0 


RUNKT+20 


STA 


0 


ST0R1 


LDA 


0 


LAM 


FMU 


0 


RUNKT+23 


F AO 


0 


ST0R1 


STA 


0 


A 1 


LAC 


0 


MU 


FMU 


0 


RUNKT+34 


STA 


0 


ST0R1 


LOA 


0 


LAM 


FMU 


0 


RUNKT+37 


FAD 


0 


ST0R1 


STA 


0 


A2 


FMU 


0 


A 1 


FMU 


0 


ADOMINWT 


FDV 


0 


HYP3 


STA 


0 


RUNKT+77 


LAC 


0 


MU 


FMU 


0 


RUNKT+50 


STA 


0 


ST0R1 


LDA 


0 


LAM 


FMU 


0 


RUNKT+53 



(-PI1) 

SET L AM 1 00T = -P I I 
(-RHO 1 ) 

SET MU1D0T = -RHO 1 

Y 

Y SQUARED 

Y SQ - X SQ 

Y SQ - 2 ( X SQ) 

(Y SQ - 2 X SQ) LAM l 

MINUS THREE 
( — 3 X ) 

(-3XY ) 

(-3XY MU1) 



SET PI 1D0T 
( — Y ) 

(- Y SQUAREO) 

(-2 Y SQUARED) 

X SQ - 2 Y SQ 
(X SQ - 2 Y SQJMUl 

MINUS THREE 
( — 3 X ) 

(-3XY ) 

I-3XY LAM 1 ) 



SET RHO 1 DOT 



( — W ) • ENTER FROM L70+5 

(-WT) 

1 — WT 

LOAD INITIAL ACCEL. 

A/l 1-WT ) 



(-MU LAM 1 ) 



LAM MU1 

A 1 * LAM MU 1 - MU LAM 1 
( -MU LAM2 ) 



LAM MU2 

A2 = LAM MU2 - MU LAM2 
A 1 TIMES A2 
(A/ 1-WTJA1A2 

SET A12D0T 

l— MU LAM3 ) 



LAM MU3 



i 



22 





FAD 


0 


STOR 1 




STA 


0 


A3 




FMU 


0 


A 1 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP 3 




STA 


0 


RUNKT+102 




LAC 


0 


MU 


L 1 70 


FMU 


0 


RUNKT+64 




STA 


0 


STOR 1 




LDA 


0 


LAM 




FMU 


0 


RUNKT+67 




FAD 


0 


STOR 1 




STA 


0 


A4 




FMU 


0 


A 1 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP3 




STA 


0 


RUNKT+ 105 




LDA 


0 


A2 




FMU 


0 


A2 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP3 




STA 


0 


RUNKT+ 1 10 




LDA 


0 


A2 


L200 


FMU 


0 


A3 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP 3 




STA 


0 


RUNKT+ 1 13 




LDA 


0 


A2 




FMU 


0 


A4 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP3 




STA 


0 


RUNKT+ 1 16 




LDA 


0 


A3 




FMU 


0 


A3 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP 3 




STA 


0 


RUNKT+ 1 2 1 




LDA 


0 


A3 




FMU 


0 


A4 


L2J0 


FMU 


0 


ADOMINWT 




FDV 


0 


HYP3 




STA 


0 


RUNKT+ 1 24 




LDA 


0 


A4 




FMU 


0 


A4 




FMU 


0 


ADOMINWT 




FDV 


0 


HYP3 




STA 


0 


RUNKT-H27 




ENI 


0 


0 




SLJ 


0 


ALPHA+2 




ENI 


0 


0 




ENI 


1 


147 




LDA 


1 


RUNKT 




STA 


1 


SAVE 




IJP 


l 


/-I 




SLJ 


4 


ALPHA 


L220 


ZRO 


0 


RUNKT 




ZRO 


0 


DERIV 




IJP 


2 


/+1 




SLJ 


0 


/+3 




ENI 


0 


0 




SLJ 


4 


ALPHA-*- 1 




ENI 


0 


0 




SLJ 


0 


/-2 




LDA 


0 


CAPXFIN 




FSB 


0 


RUNKT+4 




STA 


0 


C3 




FMU 


0 


C 3 




STA 


0 


ST0R2 




LDA 


0 


CAPYFIN 




FSB 


0 


RUNKT+7 



A3 = LAM MU3 - MU LAM3 
A 1 TIMES A3 

SET A 1 3D0T 
(-MU LAM4 ) 

LAM MU4 

A4 = LAM MU4 - MU LAM4 
A 1 TIMES A4 

SET A14D0T 
A2 SQUARED 

SET A22D0T 
A2 TIMES A3 

SET A23D0T 
A2 TIMES A4 

SET A24D0T 
A3 SQUARED 

SET A33D0T 
A3 TIMES A4 

SET A3400T 
A4 SQUARED 

SET A44D0T 

END OF DERI V PROG. 

ENTER FROM L70+7 

STORE INIT. CONDITIONS 
FOR NEXT ITERATE 

SET UP GILL ROUTINE 
PARAMETER WORD 

JUMP IF B2 NOT = ZERO 
B2=0. JUMP TO RMS COMP. 
B2 NOT = ZERO 
INTEGRATE AGAIN 

JUMP TO TEST B2 

BEGIN R.M.S. COMPUTATION 

CAPXFIN - XFIN 

CAPYFIN - YF IN 



23 



L230 



L240 



L250 



126 0 



1270 



STA 
F MU 
FAD 
STA 
LDA 
FSB 
STA 
FMU 
FAD 
STA 
LDA 
FSB 
STA 
FMU 
FAD 
SLJ 
ZRO 
FDV 
STA 
ENI 
ENI 
LDA 
STA 
INI 
IJP 
SLJ 
ZRO 
SLS 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
LDA 
STA 
STA 
LDA 
STA 
STA 
LDA 
STA 
LDA 
STA 
STA 
LDA 
STA 
ENI 
ENI 
SLJ 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 
ZRO 



0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

1 

3 

3 

1 

3 
1 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 
0 

3 

4 
0 
0 
0 
0 
0 
0 
0 
0 
0 



C4 

C4 

ST0R2 
ST0R2 
CAPUFIN 
RUNKT+ 12 
C 1 
C 1 

ST0R2 

ST0R2 

CAP VF IN 

RUNKT+15 

C2 

C2 

STOR2 

SQROOT 

0 

TWO 

RMS 

17 

55 

RUNKT+20 
BB1 1 
-3 
/-I 

MATRIX 

0 

L250-3 

0 

0 

5 

4 

B81 1 

0 

0 

0 

BB1 1+120 
RUNKT+100 
BB1 1+100 
RUNKT+ 1 03 
BB1 1+101 
RUNKT+ 106 
BB1 1+102 
RUNKT+ I 1 1 
BB1 1+103 
RUNKT+1 14 
BB1 1+104 
BB1 1+106 
RUNKT+ 1 17 
BB1 1+105 
BB1 1 + 11 1 
RUNKT+1 22 
BB1 1+107 
RUNKT+ 1 25 
BB1 1+1 10 
B81 1+112 
RUNKT+ 130 
BB1 1+113 
0 
2 

MATRIX 

L270-2 

0 

0 

4 

4 

B81 1 + 120 

4 

BB1 1+100 
3 



CAPUF IN - UFIN 



CAPVFIN - VEIN 



STORE ROOT MEAN SQUARE 

SET 4 BY 4 MATRIX 
LAMIMUIPI IRHOI = B 
(1=1,2,3,41 AT TAUFIN 

INVERT THE MATRIX 

SINGULAR MATRIX ALARM HALT 



CODE TO INVERT = 5 
L = 4 

BB 1 1 = ADDRESS OF B MATRIX 
N=0 NOT USED IN INVERT 



M = 0 
BB 1 
A 12 
SET 
A 1 3 
SET 
A 14 
SET 
A22 
SET 
A23 
SET 
SET 
A24 
SET 
SET 
A33 
SET 
A34 
SET 
SET 
A44 
SET 



NOT USED IN INVERT 
1+120 = AODRESS OF INVERSE 



TO SET 4 BY 
Dll = A 1 2 

D 1 2 = A 1 3 

D 1 3 = A 1 4 

D 2 1 = A22 

D22 = A23 
D31 = A32 

D 23 = A24 
04 1 = A42 

D32 = A33 

D 33 = A34 
D42 = A43 

D43 = A44 



3 D MATRIX 



FOR USE 6TH WORD HENCE 
MULT. B INVERSE BY D 
THIS WORD NOT USED 



4 = CODE TO MATRIC MULTIPLY 
L = 4 ROWS IN B INVERSE 
BB1 1+120 = ADDRESS OF B INVERSE 
N = 4 COLUMNS IN B INVERSE 
BB1 1+100 = ADDRESS OF D MATRIX 
M = 3 COLS. IN THE MATRIX PROD. 



24 



1300 



13 1 0 



L320 



L330 



ZRO 

REM 

REM 


0 


BB l 1+1 


LDA 


3 


BB 1 1+12 


STA 


3 


B42 


IJP 


3 


/-I 


ENI 


3 


2 


LDA 


3 


BB 1 1+7 


STA 


3 


B32 


IJP 


3 


/-I 


ENI 


3 


2 


LDA 


3 


B81 1 + 4 


STA 


3 


B22 


IJP 


3 


/-I 


LDA 


0 


RUNKT+ 1 1 


FSB 


0 


CAPUDFN 


STA 


0 


BB 1 1 


LDA 


0 


RUNKT+ 14 


FSB 


0 


CAP VDFN 


STA 


0 


B21 


LDA 


0 


RUMKT+1 2 


FSB 


0 


CAPUFIN 


STA 


0 


B31 


LDA 


0 


RUNKT+ 1 5 


FSB 


0 


CAPVFIN 


STA 


0 


B41 


SLJ 


0 


L3 1 0 


STA 


0 


BB 1 1 


ZRO 


0 


C4 


ENI 


0 


0 


SLJ 


4 


MATRIX 


SLS 


0 


L 31 0+1 


ZRO 


0 


0 


ZRO 


0 


0 


ZRO 


0 


5 


ZRO 


0 


4 


ZRO 


0 


881 1 


ZRO 


0 


0 


ZRO 


0 


0 


ZRO 


0 


0 


ZRO 


0 


BB l 1+120 


ENI 


0 


0 


SLJ 


4 


MATRIX 


ZRO 


0 


L 320-1 


ZRO 


0 


0 


ZRO 


0 


0 


ZRO 


0 


4 


ZRO 


0 


4 


ZRO 


0 


BBl 1 + 120 


ZRO 


0 


4 


ZRO 


0 


Cl 


ZRO 


0 


1 


ZRO 


0 


C 1 + 4 


ENI 


0 


0 


SLJ 


0 


/+2 


STA 


0 


BBl 1 


ZRO 


0 


Cl + 7 


ENI 


0 


0 


ENI 


0 


0 


ENI 


0 


0 


ENI 


0 


0 


ENI 


0 


0 


ENI 


0 


0 


ENI 


0 


0 


ENI 


0 


0 


IJP 


6 


/+1 


SLJ 


0 


FINAL 


SIU 


6 


/+ 1 


SLS 


2 


/+1 


ENI 


6 


0 



8B11+1 = AODRESS OF MATRIX PROO 
L273-306 SETS UP MATRIX OF SYSTEM TO B 
SOLVED FOR DELTAUFIN, VAREL , VAREM AND VA 



UDOTFIN 

UDOTFIN - CAPUDFN 


(OR 


PASS 


IF 


RENO. POINT 


FIXED 


VDOTF IN 

VDOTFIN - CAPVDFN 


(OR 


PASS 


IF 


t 

REND* POINT 


FIXED 


UFIN 











UFIN 

UFIN - CAPUFIN (OR PASS IF REND. POINT fIXED) 



VFIN 

VEIN - CAPVFIN (OR PASS IF REND. POINT FIXED) 

BYPASS PRINT OUT 
PARAMETER WORD TO PRINT 



INVERT SYSTEM MATRIX 
SINGULAR MATRIX ALARM HALT 



5 * CODE TO INVERT 
L = 4 

8B 1 1 = ADDRESS OF MATRIX 
N=0 NOT USED 

M=0 NOT USED 

BBl 1+120 = ADDRESS OF INVERSE 

MULT. INVERSE BY COL. MATRIX* 
THIS WORD NOT USED 



4 = CODE TO MATRIC MULTIPLY 

L = 4 ROWS IN INVERSE 

BB1 1+120 = ADDRESS OF INVERSE 

N * 4 COLS. IN INVERSE 

Cl * ADDRESS OF COLUMN MATRIX 

M=1 COL. IN THE COL. MATRIX 

C 1+4= ADD. OF COL. MATRIX RESULT 

BYPASS PRINT OUT 
PARAMETER WORD TO PRINT 



JUMP IF B6 NOT = ZERO 

B6 = 0. JUMP TO FINAL PROG. 

STOP SWITCH 
RESET B6 



- 25 - 





6NI 


1 


3 






LOA 


I 


C 1 + 4 






STA 


1 


VARTAU 


STORE DELTAUFIN, VAREL, VAREM.VAREN 




IJP 


1 


/-I 






FAD 


0 


TAUFIN 


NEW TAUFIN FOR NEXT ITERATE 




STA 


0 


SAVE+143 


STORE NEW TAUFIN (NOTE L30+3 ) 




ENI 


0 


0 


OR SLJ 0 L350+5 IF REND. POINT FI 


L340 


FMU 


0 


W 






FAD 


0 


S 






STA 


0 


SPWT 






SLJ 


4 


TRIG+70 






STA 


0 


COSSPWT 






FMU 


0 


R 






STA 


0 


SAVE+136 


NEW CAPYFIN 




LDA 


0 


SPWT 






SLJ 


4 


TRIG 






ZRO 


0 


0 






STA 


0 


S IN SPWT 






FMU 


0 


R 






STA 


0 


SAVE+ 135 


NEW CAPXFIN 




LDA 


0 


COSSPWT 






FMU 


0 


V 






STA 


0 


SAVE+ 137 


NEW CAPUFIN 


1350 


LAC 


0 


SINSPWT 






FMU 


0 


V 






STA 


0 


SAVE+140 


NEW CAPVFIN 




LAC 


0 


COSSPWT 






FMU 


0 


VSQDR 






STA 


0 


SAVE+142 


NEW CAPVDFN 




LAC 


0 


SINSPWT 






FMU 


0 


VSQDR 






STA 


0 


SAVE+141 


NEW CAPUDFN 




ENI 


0 


0 






ENI 


0 


0 






ENI 


1 


2 






LDA 


1 


C 1 + 5 






FAD 


1 


EL 






STA 


1 


EL 


STORE NEW EL, EM AND EN 




IJP 


1 


/-I 




L360 


SLJ 


4 


700C7 


PRINT 




ZRO 


0 


0 






SAU 


0 


MDOT 


OECIMAL DUMP WITH PANEL 




ZRO 


0 


S 


i 




EXF 


0 


70 


CLEAR ARITH. ERRORS 




SLJ 


4 


70007 


PRINT 




STA 


0 


CAPXFIN 


DEC. DUMP 




ZRO 


0 


LAM 






SLJ 


0 


RECUR 


TO NEXT ITERATION 




ZRO 


0 


0 






BSS 




3 




POLAR 


BSS 




200 


B 1 = 33 B 6=6637 0 P=460 


EL 1 


BSS 




1 




EMI 


BSS 




1 




ENT 


BSS 




1 




T 1 


BSS 




1 




EL2 


BSS 




1 




EM2 


BSS 




1 




EN2 


BSS 




1 




T2 


BSS 




1 




EL3 


BSS 




1 




EM3 


BSS 




1 




EN3 


BSS 




1 




T3 


BSS 




1 




EL4 


BSS 




1 


• 


EM4 


BSS 




1 




EN4 


BSS 




1 




T4 


BSS 




1 




EL5 


BSS 




1 




EM5 


BSS 




1 




EN5 


BSS 




1 





26 



T5 


BSS 




1 






BSS 




26 


STORE FOR MATRIX PRODUCT 


FINAL 


ENI 


1 


147 


COMPUTES POINTS ON FINAL TRAJ. 




LIL 


2 


FIXIPLON 






LDA 


1 


SAVE 






STA 


1 


RUNKT 


RESET INITIAL CONDITIONS 




IJP 


1 


/-I 






SLJ 


4 


ALPHA 


SET UP GILL ROUTINE 




2R0 


0 


RUNKT 


PARAMETER WORD 




ZRO 


0 


DERIV 






ENI 


0 


0 






SLJ 


4 


/ +4 


RET. JUMP TO DUMP SUBROUTINE 




IJP 


2 


/+ 1 


JUMP IF B2 NOT = ZERO 




SLJ 


0 


L70C+1 


B2=0. JUMP TO REMAINDER OF FINAL PROG 


L650 


ENI 


P 


0 






SLJ 


4 


ALPHA+1 


INTEGRATE AGAIN 




ENI 


Q 


0 


4 




SLJ 


0 


/-3 






SLJ 


0 


0 


ENTRY TO DUMP SUBROUTINE ' 




LDA 


0 


RUNKT +2 


TAU 




ENI 


0 


0 






ENI 


0 


0 






STA 


0 


DUMP 


(DUMP) = TAU 




LDA 


G 


RUNKT+4 


X 




ENI 


0 


0 






ENI 


0 


0 






STA 


0 


DUM P+4 


(DUMP+4) = X 




LDA 


0 


RUNKT +7 


Y 




ENI 


0 


0 






ENI 


0 


0 




L660 


STA 


0 


DUM P+5 


(0UMP+5I = Y 




LDA 


0 


RUNKT+12 


U 




ENI 


0 


0 






ENI 


0 


0 






STA 


0 


DUM P+6 


(DUMP+6) = U 




LDA 


0 


RUNKT+15 


V 




ENI 


0 


0 






ENI 


0 


0 






STA 


0 


DUMP+7 


(DUMP+7) = V 




LDA 


0 


RUNKT+1 1 


UDOT 




ENI 


0 


0 






ENI 


0 


0 






STA 


0 


DUMP+10 


(DUMP+10) = UDOT 




LDA 


0 


RUNKT+ 14 


VDOT 




ENI 


0 


0 






ENI 


0 


0 




L670 


STA 


0 


DUMP+1 1 


(DUMP+11) = VDOT 




LDA 


0 


RUNKT+12 


U 




LDQ 


0 


RUNKT+15 


V 




SLJ 


4 


POLAR+ 1 30 






ENI 


0 


0 






ENI 


0 


0 






STQ 


0 


DUMP+12 


(DUMP+12) = ARCTAN ( V/U ) 




LDA 


0 


LAM 






LDQ 


0 


MU 






SLJ 


4 


POLAR+ 130 






ENI 


0 


0 






ENI 


0 


0 






STQ 


0 


DUMP+13 


(DUMP+13) = ARCTAN (MU/ LAM) 




SLJ 


4 


70007 


PRINT 




STA 


0 


DUMP 


PARAMETER WORD 




ZRO 


0 


DUMP+13 




L700 


SLJ 


0 


L 650+2 


TO EXIT OF DUMP S.R. 




ZRO 


0 


0 




LDA 


0 


C3 


C3 = CAPXFIN - XF IN 




ENI 


0 


0 






ENI 


0 


0 






STA 


0 


OUMP+14 


(OUMP+14) = CAPXFIN - XF I N 




LDA 


0 


C4 


C 4 = CAPYFIN - YF IN 




ENI 


0 


0 



27 





ENI 


0 


0 




STA 


0 


DUMP+ IS 




LDA 


0 


C 1 




ENI 


0 


0 




ENI 


0 


0 




STA 


0 


DUMP+16 




LDA 


0 


C2 




ENI 


0 


0 


L710 


ENI 


0 


0 




STA 


0 


DUMP+17 




SLJ 


4 


700C7 




ZRO 


0 


0 




STA 


0 


DUMP+14 




ZRO 


0 


DUMP+17 




SLJ 


0 


L3570+5 




ZRO 


0 


0 


STOR1 


BSS 




1 


ST0R2 


BSS 




1 


STOR3 


BSS 




1 


ST0R4 


BSS 




1 


L720 


DEC 




1. 


DELX1 


BSS 




1 


DELY 1 


BSS 




1 


DELU1 


BSS 




1 


DEL V 1 


BSS 




1 




DEC 




1 . 


DELX2 


BSS 




1 


DEL Y2 


BSS 




1 


DELU2 


BSS 




1 


DELV2 


BSS 




1 




DEC 




1 . 


DELX3 


BSS 




1 


DELY3 


BSS 




1 


DELU3 


BSS 




1 


DELV3 


BSS 




1 




DEC 




1. 


DELX4 


BSS 




1 


DELY4 


BSS 




1 


DELU4 


BSS 




1 


DELV4 


BSS 




1 




DEC 




1. 


DELX5 


BSS 




1 


DEL Y5 


BSS 




1 


DELU5 


BSS 




1 


SAVE 


BSS 




150 


RUNKT 


BSS 




135 


CAPXFIN 


BSS 




1 


CAPYFIN 


BSS 




1 


CAPUFIN 


BSS 




1 


CAPVFIN 


BSS 




1 


CAPUDFN 


BSS 




1 


CAPVDFN 


BSS 




1 


TAUFIN 


BSS 




1 


LAM 


BSS 




1 


HYP 


BSS 




1 


HYP3 


BSS 




1 


MU 


BSS 




1 


HYP 2 


BSS 




1 


UNITY 


DEC 




T 

1 • 


W 


BSS 




1 


DELTAU 


BSS 




1 


SPWT 


BSS 




1 


COSSPWT 


BSS 




1 


SINSPWT 


BSS 




1 


R32 


BSS 




1 


R52 


BSS 




1 


IMTSIGN 

FRACSIGN 


BSS 




1 


BSS 




1 


FLOIPLON 


BSS 




1 


FIXIPLON 


BSS 




1 



( DUMP + 15) = CAPYFIN - 
Cl = CAPUFIN - UF IN 



{ DUMP + 16) = CAPUFIN - 
C2 «* CAPVFIN - VF IN 

(DUMP+17) = CAPVFIN - 
PRINT 

PARAMETER WORD 

JUMP TO PRINT B FINAL 



ALSO DELV5 STORE 



- 28 - 



YFIN 
UF I N 
VF I N * 



THREDR2 


OCT 




2002417 


TREHALF 


DEC 




1.5 


DELTFLSG 


BSS 




1 


NT 


BSS 




l 


THETA 


BSS 




1 


VSTART 


DEC 




.65 


DELTHETA 


BSS 




1 


DELXFIN 


BSS 




1 


DELYFIN 


BSS 




1 


DELUFIN 


BSS 




1 


DEL VF IN 


BSS 




l 


MOOT 


DEC 




.0025 


A 


BSS 




1 


QQ 


BSS 




1 


C 


DEC 




10000. 


EL 


BSS 




1 


EM 


BSS 




1 


EN 


BSS 




1 


S 


BSS 




1 


A1 


BSS 




1 


A2 


BSS 




1 


A3 


BSS 




1 


A4 


BSS 




1 


VARTAU 


BSS 




l 


VAREL 


BSS 




1 


VAREM 


BSS 




1 


VAREN 


BSS 




1 


ALPHA 


BSS 




134 




BSS 




14 


DUMP 


BSS 




20 


BB! 1 


BSS 




1 


B 1 2 


BSS 




1 


B 1 3 


BSS 




1 


B 1 4 


BSS 




1 


B21 


BSS 




1 


B22 


BSS 




1 


B23 


BSS 




1 


B24 


BSS 




1 


B31 


BSS 




1 


832 


BSS 




1 


B33 


BSS 




1 


B34 


BSS 




1 


B41 


BSS 




1 


B42 


BSS 




1 


B4 3 


BSS 




1 


B44 


BSS 




1 




BSS 




120 


Cl 


BSS 




1 


C2 


BSS 




1 


C3 


BSS 




1 


C4 


BSS 




1 




BSS 




4 


LOG 


BSS 




100 




REM 








REM 






L2000 


BSS 




1440 




REM 






L3440 


ENA 


0 


0 




ENI 


1 


17 




STA 


1 


BB 1 1 




IJP 


1 


/ 




LDA 


0 


TAUFIN 




SLJ 


4 


TRIG 




STA 


0 


BB 1 1 




LDA 


0 


TAUFIN 




SLJ 


4 


TRIG+70 




ZRO 


0 


0 




STA 


0 


B21 




LDA 


0 


TAUFIN 




FMU 


0 


R00T2 



6663160 THREE /R00T2 



NON-DIMEN. INITIAL SPEED 



FRACTIONAL MASS LOSS PER SEC. 



ROCKET NOZZLE VEL. IN FT. /SEC, 



INI T. TARGET SECTOR ANGLE 



B 1 *30 B6=67340 P=460 ' 



5 TIMES L SQUARED CELLS 



STORE FOR MATRIX PRODUCT 
B 1 = 12 B6=67700 P=460 

LOG S.R. ALARM EXIT TO K110 
(67703) = 750 65110 000 00000 
USNPGS GEN. DUMP Bl=32 P=2 1 

AID TO LINEARIZED TRAJECTORY PROG< 



CLEAR BB I 1 MATRIX AREA 



(BB11 ) ■ SINT- 



(B21) = COST 



29 



SINHTR2 



L3450 



L 3460 



L3470 



L 3500 



L3510 



SLJ 


4 


EXP+7 1 




STA 


0 


832 


( B32 ) = SINHTR2 


LDA 


0 


TAUFIN 




FMU 


0 


R00T2 




SLJ 


4 


EXP+107 




STA 


0 


B42 


( B42) = C0SHTR2 


LDA 


0 


SIN 


SINE IMPULSE INTEGRAL 


FDV 


0 


C 




STA 


0 


813 




LDA 


0 


COS 


COSINE IMPULSE INTEGRAL 


FDV 


0 


C 




STA 


0 


823 




LDA 


0 


SINH 


SINH IMPULSE INTEGRAL 


FDV 


0 


C 




STA 


0 


833 




LDA 


0 


COSH 


COSH IMPULSE INTEGRAL 


FDV 


0 


C 




STA 


0 


R43 




LAC 


0 


CAP YF IN 




STA 


0 


B 14 




LAC 


0 


CAPVFIN 




STA 


0 


824 




SLJ 


0 


L357C+3 


JUNMP TC PATCH 


STA 


0 


B34 




LDA 


0 


CAPUFIN 




STA 


0 


844 




LAC 


0 


DEL XF IN 




STA 


0 


Cl 


Cl « ADD. OF MATRIX COL. 


LAC 


0 


DELUFIN 




STA 


0 


C2 




LAC 


0 


DELYFIN 




FMU 


0 


R00T2 




STA 


0 


C3 




LAC 


0 


DELVFIN 




STA 


0 


C4 




ENA 


0 


/+2 




SAL 


0 


L320+4 


CREATE LINEAR EQU. ROUTINE EXIT 


SLJ 


0 


L310 


JUMP TO LINEAR EQU. ROUTINE 


ZRO 


0 


0 




ENA 


0 


L 320+6 




SAL 


0 


L320+4 


REBUILD RENDEVOUS PROGRAM 


LOA 


0 


C 1 + 4 


DEL VI 


FAD 


0 


V 1 




STA 


0 


ST0R1 


(ST0R1 ) * NEW VI 


FMU 


0 


ST0R1 




STA 


0 


ST0R2 


( S TOR 2 ) = NEW VI SQUARED 


LDA 


0 


C2 + 4 


DEL V 2 


FAD 


0 


V2 




STA 


0 


ST0R3 


(ST0R3) = NEW V2 


FMU 


0 


ST0R3 




FAD 


0 


ST0R2 


NEW VI SQ. + NEW V2 SQ. 


SLJ 


4 


SQROOT 




ZRO 


0 


0 




STA 


0 


VST ART 1 


NEW VSTART 


LDA 


0 


ST0R1 


NEW VI 


FDV 


0 


VSTART1 


COSINE NEW THETA 


SLJ 


4 


POLAR 


TO INVERSE COSINE S.R. 


STA 


0 


THETA 1 


NEW THETA 


LDA 


0 


C3+4 


DEL C 


FAD 


0 


C 


NEW C 


STA 


0 


EJECT 


(EJECT) = NEW C 


LDA 


0 


C4+4 


DIFFERENTIAL OF B 


FAD 


0 


B 


NEW B 


STA 


0 


SECTOR 


( SECTOR) = NEW B 


SLJ 


4 


70007 


PRINT RESULTS DECIMAL 


STA 


0 


VST ART 1 


PARAMETER WORD 


ZRO 


0 


SECTOR 




SLJ 


4 


70007 


PRINT RESULTS OCTAL 


ZRO 


0 


0 




ZRO 


0 


VST ART 1 


PARAMETER WORD 



- 30 





ZRO 


0 


SECTOR 






ENI 


r, 

V* 


0 






ENI 


0 


0 






LOA 


0 


VSTART 1 


RESUME COMP. OF INIT. CONDITIOI 




STA 


0 


V START 


SET NEW VSTART 




LOA 


0 


THE TA 1 






STA 


0 


THETA 


SET NEW THETA 




LDA 


0 


EJECT 






STA 


0 


C 


SET NEW C 




LDA 


0 


SECTOR 






STA 


0 


B 


SET NEW 8 




SLJ 


4 


TRIG+70 






ZRO 


0 


0 




L3520 


STA 


0 


COSB 






FMU 


0 


R 






STA 


0 


CAPYFIN 


TARGET COOR. AT RENDEVOUS 




LDA 


0 


B 






SLJ 


4 


TRIG 






ZRO 


0 


0 






STA 


0 


SINB 






FMU 


0 


R 






STA 


0 


CAPXFIN 


TARGET COOR. AT RENDEVOUS 




LDA 


0 


C0S8 






FMU 


0 


V 






STA 


0 


CAPUFIN 


TARGET VEL. COMP. AT RENDEVOUS 




LAC 


0 


SINB 






FMU 


0 


V 






STA 


0 


CAPVFIN 


TARGET VEL. COMP. AT RENDEVOUS 




LAC 


0 


COSB 




L3530 


FMU 


0 


VSQDR 






STA 


0 


CAP VDFN 


TARGET ACC. COMP. AT RENDEVOUS 




LAC 


0 


SINB 






FMU 


0 


VSQDR 






STA 


0 


CAPUDFN 


TARGET ACC. COMP. AT RENDEVOUS 




LDA 


0 


THETA 






SLJ 


4 


TRIG+70 


COSINE THETA 




ZRO 


0 


0 






FMU 


0 


VSTART 






STA 


0 


V 1 


INITIAL HORIZONTAL VEL. 




LDA 


0 


THETA 






SLJ 


4 


TRIG 


SINE THETA 




FMU 


0 


VSTART 






STA 


0 


V2 


INITIAL VERTICAL VEL. 




LDA 


0 


C 


ROCKET NOZZLE VEL. 




FSB 


0 


TEN GRAND 


NOZZLE VEL. DIFFERENCE 


L3540 


A JP 


2 


/+1 


JUMP IF DIFF. IS POS. 




SCM 


0 


MASK 


ABSO. VALUE OF DIFF. 




THS 


0 


UNITY 


EXIT IF ONE GREATER THAN DIFF. 




SLJ 


0 


SEC TOR+2 


NO. JUMP TO ADJUST MOOT AND C 




LDA 


0 


C 


YES. LOAD ROCKET NOZZLE VEL. 




FMU 


0 


MOOT 






FDV 


0 


GACCEL 






STA 


0 


A 


NON-DIMEN. ACCELERATION 




LDQ 


0 


ELS 






SLJ 


0 


K10C+2 


RETURN TO LIN. TRAJ. PROG. 


VSTART 1 


BSS 




1 


TEMP. STORE FOR NEW VSTART 


THETA 1 


BSS 




1 


TEMP. STORE FOR NEW THETA 


EJECT 


BSS 




I 


TEMP. STORE FOR NEW C 


SECTOR 


BSS 




1 


TEMP. STORE FOR NEW B 




SLS 


0 


SECTOR+1 


HALT IF 1-WT IS NEGATIVE 




ZRO 


0 


0 






LAC 


0 


TAUFIN 


(-T) 




FMU 


0 


OMEGA 


(-WT) . 




FAD 


0 


UNITY 


( 1-WT ) 




AJP 


3 


SECTOR+1 


HALT IF 1-WT IS NEGATIVE 




SLJ 


4 


LOG+66 


LOG(l-WT) TO BASE E 




ZRO 


0 


0 






FMU 


0 


C 






FDV 


0 


TENGRAND 






SLJ 


4 


EXP 





31 





ZRO 


0 


0 




SCM 


0 


MASK 




FAD 


0 


UNITY 


L3560 


FOV 


0 


TAUFIN 




STA 


0 


OMEGA 




LDA 


0 


GACCEL 




FDV 


0 


REARTH 




SLJ 


4 


SQROOT 




ZRO 


0 


0 




F MU 


0 


OMEGA 




STA 


0 


MDOT 




LDA 


0 


TENGRAND 




STA 


0 


C 




ENI 


0 


0 




ENI 


0 


0 




ENI 


0 


0 




ENI 


0 


0 




ENI 


0 


0 




ENI 


0 


0 


L3570 


ENI 


0 


0 




LDA 


0 


MDOT 




FMU 


0 


C 




FDV 


0 


GACCEL 




STA 


0 


A 




SLJ 


0 


K70 




LDA 


0 


CAPXF IN 




FMU 


0 


ROOT 2 




SLJ 


0 


L34Z0+2 




ZRO 


0 


0 




SLJ 


4 


70007 




ZRO 


0 


0 




STA 


0 


UNITY 




ZRO 


0 


SPWT 




SLS 


0 


L 3570+7 




ZRO 


0 


0 




END 







NEW OMEGA 
STORE NEW OMEGA 



SHULER ANGULAR FREQ. 

NEW MOOT 
STORE NEW MDOT 
LOAO NEW C 
STORE NEW C 



STORE NEW A 

JUMP TO LIN. TRAJ. PROG. 
PATCH FCR L3U60+1 L.I. 



ENTRY FROM L710+3 

PRINT B FINAL DECIMAL 

PARAMETER WORD 

FINAL B = SPWT 

END OF NON-LINEAR PROGRAM 
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•U6 Bleick 
no. 34 Orbital transfer in 

minimum time. 
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Bleick 

34 Orbital transfer in 

minimum time. 
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